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SUMMARY 


Energy levels and electron densities are calculated in different oxidation states of diphosphopy- 
ridine nucleotide and flavine-adenine dinucleotide, both of which take part in the oxidative 
processes in the living cell. Only those parts of the molecules where oxidation and reduction take 
place are considered in the calculation. The z-electrons and some lone pairs in these parts are 
considered explicitly. Molecular orbitals (MO) for these “‘z-electrons’’ are formed by linear 
combination of atomic orbitals (LCAO). The antisymmetrized products (AP) of these orbitals 
are used. Some of the integrals in the electronic wave equation are determined semi-empirically 
according to a treatment by Pariser and Parr [1] and by Pople [2]. The self-consistent field method 
(SCF) is used for solving the equations. In the discussion it is shown that the results seem to give 
considerable support to the hypothesis of the present writer [3] concerning the oxidative phos- 
phorylation at the diphosphopyridine nucleotide and flavine molecules in the living cell. 


I. Introduction 


In the present paper energy levels and electron distributions have been calculated 
in different oxidation states of diphosphopyridine nucleotide and flavine-adenine 
dinucleotide which take part in oxidative processes in the living cell. Most substances 
in the living cell are not oxidized directly by oxygen. Instead there is a chain of 
carriers which transfer electrons or hydrogen atoms from the substance to oxygen. 
These carriers alternate between the oxidized and the reduced states, thus being 
reduced by the nearest neighbour on the substrate side and oxidized by the nearest 
neighbour on the oxygen side. Diphosphopyridine nucleotide and flavine-adenine 
dinucleotide are known to be neighbours in such a chain, diphosphopyridine nucleo- 
tide being oxidized by the flavine. They are very large molecules. However, only the 


Abbreviations used: DPN+ and DPNH diphosphopyridine nucleotide in oxidized and reduced 
form; DPNH+ the free radical obtained if an electron is removed from DPNH; FAD and FADH, 
flavine-adenine dinucleotide in oxidized and reduced form and FAD~ the free radical obtained 


if an electron is transferred to FAD. 
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part of the molecule which changes oxidation state has been included in the present 
considerations, i.e. the pyridine ring of the diphosphopyridine nucleotide and the 
isoalloxazine ring of the flavine. The structure of the actual parts of the compounds 
is given in Fig. 1, the rest of the molecules being denoted R. The abreviations DPN*, 
DPNH, FAD and FADH, given in the figure are used in the following. Since it has 
been suggested in [3] that an electron is transferred from DPNH to FAD at an 
intermediate state in the oxidation in the cell, the free radicals DPNH+ and FAD— 
have been treated also, although it is most probable that there is a charge transfer 
complex between DPNH and FAD built in the cell and no DPNH* and FAD~ free 
radicals. 

Pullman and Pullman [4] have previously calculated the electron distribution and 
energy levels in different oxidation states of the diphosphopyridine nucleotide and 
of flavine mononucleotide, the latter having the same flavine part as FAD. They have, 
however, used other methods and parameter values than the present author. Their 
calculations include more electrons than the present, but on the other hand they have 
used a more approximate method, viz. the Hiickel approximation. It will therefore 
be of interest to compare the results of the different calculations. (See Section IV.) 

In Section II a description of the calculation is given. Thus IT a) includes an outline 
of the method used and states which atoms are included in the calculation and which 
electrons are considered explicitly. In II b) the wave equation for the z-electrons is 
given and II c) includes parameter values chosen according to experimental data and 
a description of the evaluation of some integrals. In II d) the method for the solution 
of the wave equation is given and in II e) and f) equalities for the calculations of 
electron densities, bond orders and total energies are given. The results are given in 
Section III and discussed in Section IV, where the hypothesis in [3] is discussed with 
respect to the present results. 


II. Calculation 


a) An outline of the method used 


The method used for the calculation is a semi-empirical LCAO-MO—AP-SCF 
method. This implies that molecular orbitals (MO) have been built by linear combina- 
tion of atomic orbitals (LCAO) and that the ordinary method with antisymmetrized 
products (AP) of MO:s has been used according to Roothaan [5a]. For the free 
radicals with unparied electrons the extended calculation by Roothaan [56] has 
been followed. Instead of calculating all the integrals a semi-empirical method given 
by Pariser and Parr [1] and by Pople [2] has been used. However, no configuration 
interaction has been taken into consideration as Pariser and Parr have done. For 
solution of the equations the self-consistent field (SCF) treatment has been used. 
More details are given in the following. 

As is mentioned in Section I, the calculation includes only parts of the large 
molecules, i.e. only the atoms numbered in Fig. 1 are considered. 

The usual z-approximation has been used, i.e. the interaction of all the z-electrons 
and lone pairs at some of the N and O atoms have been considered explicitly, the 
effect of the o-electrons being included in the term of the Hamiltonian denoted H°°"°. 
Thus for DPN* 6 electrons have been considered explicitly (i.e. the z-electrons), for 
DPNH 8 electrons (i.e. the z-electrons plus the lone pair at the N atom plus 2 elec- 
trons in the C-—H, bond, the H, group being treated as a unity [6]), for FAD 10 
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Fig. 1. The numbers give the interatomic distances, in A, used in the calculations. They are 
estimated from distances in similar bonds, given in the literature. 


electrons (i.e. 8 z-electrons plus the lone pair at the N atom denoted (7) in Fig. 1) and 
for FADH, 12 electrons (i.e. 6 z-electrons plus the lone pairs at N (7), N(8) and O(9) 
in Fig. 1). The free radical DPNH* has been calculated as DPNH, where one elec- 
tron from the occupied orbital of highest energy has been removed, and FAD-, as 
FAD in which one electron has been transferred to the empty z-orbital of lowest 
energy from a reducing substance. In these cases, therefore, 7 electrons and 11 elec- 
trons respectively have been considered explicitly. All the electrons considered 
explicitly are denoted z-electrons. 


b) The wave equation for the z-electrons 


To obtain MO:s for the z-electrons atomic orbitals v, (ju), where the index p indi- 
cates the atom to which the orbital belongs and y indicates the electron considered, 
have been combined to orthonormal molecular space orbitals, p;(/) 


Qi (Me) = Cin Xr (MH) + CinXe(M) H+ + Cis Hs (4), (1) 
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where 1, 2, ... s includes all the atoms considered in the calculation. These @; (1) have 
finally been antisymmetrized in the usual manner, giving the wave function 


ae 
ADAM ta, once) ae > P gi (1) 9g; (2), .-. @v (n) x spin functions, (2) 
n! 


where n indicates the number of electrons considered explicitly. 
The total Hamiltonian operator H for the electrons has been expressed as 


H= Dy ee 


lt — - Vv=ye 


, (3) 


where the summation over ju and y has to be carried out over all the electrons consid- 
ered explicitly. The expression H°°'* (1) is the kinetic energy operator of the electron 
u, T'(u), plus the potential energy operator of the core, U°"° (wz), U%'* (iu) being the 
field from all the nuclei plus the inner electrons plus the o-electrons. 


Hore (4) = T (u)+ pyoore (n)= FEUD Els 2 (Wp (4), (4) 


where the summation over s has to be carried out over all the atoms considered. 
In the cases where all occupied space orbitals contain two electrons with opposite 
spin the problem is to solve the equation 


Le (te) 2 2 Ji (u) — 2 K;(u)) gi (4) = &: gi (uw), (5) 
where 
Ji (UW) p(w) = lies (v) — woyde,| p (1) (6) 
and 


2 


K;(u) p (u) = i Vi (”) . vorde,| i (Ut). (7) 
uy 

The summation over 7 in (5) has to be carried out over all the occupied space orbitals. 

—é, is the ionization energy for removal of an electron in qj. 

For the calculation of the electron distribution in DPNH+ and FAD-, where 
one electron is unpaired, the method of Roothaan [5b) has been used. In applying 
his method to the present case it is found that the problem is to solve the following 
equations 


core < a 1 ie 
E (u) +2 7 J(u) — 2 K;(u) Piel) 5 Ke (u)+ UM. v| P(e) =ExG~u(u) (Sa) 
for the paired electrons and 


J 


[He (14) + 2D; (u) — > Kj (uw) + 23M; (1) Pe (Mt) = Ee Pe (Mr) (8b) 


for the unpaired electron. J; and K; are defined as in (6) and (7) and 
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M;(u) 9 (u) = 


9 


f Vi (A) pe (v) : P (?) Pe (A)dvjd 7 Pi (uw) + 


+ f pF (A) @ (A) dey J ee (v) — Qj (¥) dv, Pe wn}. (9) 
uv 


The summation over 7 in (8) has to be carried out over all the doubly occupied space 
orbitals. The index e denotes the singly occupied orbital and J,, K, and M, are 
defined in accordance with (6), (7) and (9). In the cases calculated in this paper —€ 
is approximately equal to the ionization energy, although, according to Roothaan 
[56], this is not true in the general case, when there are unpaired electrons. 


c) Parameter values and evaluation of some integrals 


To solve (5) and (8) it is necessary to know the interatomic distances in the mole- 
cules. The interatomic distances between nearest neighbours used in the calculations 
are given in Fig. 1. It should be pointed out that there seems not to be any experi- 
mentally determined distances given in the literature, as far as the author is aware, 
so that the distances given in the figure are estimated by comparison with distances 
in similar bonds. For the free radicals DPNH+ and FAD~- the same interatomic 
distances as in DPNH and FAD have been used. 

The first step in the calculation has been to determine the values of the following 
integrals. 


ap = | Xp (mu) A? (u) ky (Mw) dey (10) 
Boa=J Xp (u) A" () Lq(u) dr (11) 
and 
2 
(pr|qs)=J Xp (u) Za (r) = Ar (u) Xs (v) dv, d v,. (12) 


pv 


Following Pariser and Parr [1], (10) and (11) are assumed to be parameters, to be 
determined semi-empirically. Furthermore, the calculation has been simplified by 
assuming tentatively that #,, =0, when p and q are non-neighbours, and that (pr | qs) 
=0, when p+r and q +s, ie. only integrals of the type (pp/qq) are taken into 
consideration. The last assumption implies a formal neglect of differential overlap, 
ie. yxy,=0 for p+r. Furthermore, in a, the terms f75(u) UP" (u) Z(t) avy (see 
(10) and (4)), when p +q, have been calculated as 


$5 (ue) US"? (12) hp (i) by = f Lp (wu) UG (1) Xv (1) dy, — Na (pP|d9): 


where N, = 1 or 2, according to the number of electrons contributed to the z-system 
by the atom q. The first integral in the right-hand member will always be small. 
Following Pople [2] we have neglected this integral in the present computations. Thus 


tp = —I, ~ 2 (ppl 99), (13) 
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where — I, is the ionization potential of the valence state of the free atom. The summa- 
tion over g has been carried out over all the electrons considered explicitly. 

Following Pariser and Parr [1] and Anno [7] the ionization potentials and the 
f-values have been taken from empirical data. The values used are given in Table 1. 
The J, and Jy have been taken from Anno [7] and Jy, chosen = 26.2 eV. I, has been 
taken from Skinner and Pritchard [8] and Iy, from [9]. The 6-values are dependent 
on the distance between the atoms. For Boo the values have been taken from Pariser 
and Parr [1], who have given fgc¢-values for two distances and a formula for deter- 
mining values for other distances. For cy the value — 2.74 eV for the distance 1.35 A 
has been used, in accordance with Anno’s suggestions [7]. To get /oy-values for 
other distances, it has been assumed that the curve giving the dependence of Poy on 
the interatomic distance is parallel to the corresponding curve for foc. boo for the 
distance 1.23 A has been taken from Sidman [10]. For the distance 1.40 A it has been 
supposed that Boo (7 = 1.40 A) = Box (v7 = 1.40 A). The value of fox, has been obtained 
by calculating the integrals in 


Bav = Apa” — (Spq/2) (Hpn* + Haq"). (14) 


See footnote in [la]. For boy. the value —3.6 eV has been chosen. 

The expressions (pp/qq) have been determined according to formulas given by 
Pariser and Parr [1]. For interatomic distances r> 2.80 A they have given the 
following formula obtained from the uniformly charged sphere model: 


(pp|qq)=7.1975/r{[1 + (1/27)? (Rp — Be)" + [1+ (/2r) (Rp + Ra)? }eV, (15) 
R, and R, being diameters of nonconducting spheres. They are to be determined as 
R, = 4.597 /Z, (in A), (16) 

Z, being Slater’s effective nuclear charge. To get a corresponding value on Ry, the 
present author has determined the diameter Ry, by comparing the ball values and 


the exact values of integrals at 7 =0 (compare Anno [7]). All the R, values used in 
the present paper are given in Table 2. 


Table I. Values of ionization energies, J, and of resonance integrals, f, (11), used 
for calculations. 


Ig Ty Iy+ Ty, Io 
IE a oe ee 
eV 11.54 14.00 26.2 15.6 We/eealk 
evatte 0.4242 0.5146 0.962 0.573 0.6326 
_ 
B 
Interatomic sa we Oe Poo ie 
distances, A 1.34 1.39 1.47 1.33 1.36 1.38 1.42 PR} 1.40 1.09 
eee 
eV = sill} 239) = Iai — 3.05 —2.60 -—2.37 —1.95 3.0 
‘ : : Ws : 4 = 8h —1.95 — 2.64 
a.u. — 0.1162 —0.0879 —0.0555 -—0.1121 —0.0956 —0.0871 —0.0717 —0.110 -—0.0717 —0.09¢ 


ee ee 
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Table 2. Diameters used in (15). 


C N N+ @) H, 


R 1.41 1.18 1.08 1.01 1.8 


For the determination of (pp/qq) for r< 2.80 A Pariser and Parr [1] have given 
the following formula for extrapolation: 


(pp|qq) = 1/2{(pp| pp) + (qq\aq)] — ar — br?, (17) 


where (pp|pp) and (gq|qq) are to be obtained from experimentally determined ioni- 
zation energies, I, of valence states of the free atom and electron affinities, A, as 


(ss|ss) =I, —A,. (18) 


The parameters a and } are to be determined from values of (pp | qq) for large r values, 
e.g. for r = 2.80 A and r = 3.70 A, calculated according to (15). 

The parameter values used in (17) in the present paper are given in Table 3 for 
different atoms p and q. 


Table 3. Formulas for determination of (pp|qq) in eV, according to (17), when the 
interatomic distance r<2.80 A. 


CC) =11.08 — 2.9226 r+ 0.26203 r? (According to Anno [7]) 
| NN) =(11.08 + 13.46)/2 — 3.6807 r+ 0.3747 r? 
| N+N+) =(11.08 + 17.21)/2 —4.841 r+ 0.5518 r? (17.21 according to Anno [7]) 
Clsk = (11.08 + 15.62)/2 — 4.4708 r+ 0.5054 r? 
CC| OO) = 12.757 — 3.966 7+ 0.4169 72 (According to Sidman [10]) 
NN | NN) =13.46 — 4.3907 r + 0.4813 r? 

| ) =(13.46 + 14.67)/2 —4.741 r+ 0.5323 r? (14.67 according to Sidman [10]) 


d) The solution of the electronic wave equation 


To start the SCF calculation a zero-order approximation of the ¢ coefficients of 
(1) are needed for the insertion in J, K and M in (5) and (8). To get such values a 
calculation in the Hiickel approximation has been carried out, i.e. all the electrons 
are included in H®”®, so that all the e?/r,,, terms in (3) disappear. In each cycle of the 
SCF calculation the solving of the system of equations corresponding to (5) or (8) to 
get the ¢ and c values has been made on the Swedish electronic computer BESK. 
The calculation has been continued until the difference between the ¢-values in the 
last two cycles are at most 0.002 a.u., ensuring at Jeast two and in most cases three 
decimals in the electron densities calculated from the ¢ values, 


e) Determination of electron densities, charges on the atoms and bond orders 


The final ¢ values in (1) have been used to calculate the electron densities at the 
different atoms according to Coulson [11]. His definition of “charge” is identical 


103 


BARBRO GRABE, Diphosphopyridine nucleotide 


with Mulliken’s definiton [12] of gross atomic population V(r) on the atom 7 when 
the overlap integrals are neglected: 


N (r) => Nich. (19) 
j 


The summation should be carried out over the occupied space orbitals and NV; should 
be 2 for doubly occupied orbitals and 1 for singly occupied orbitals. 
The charge Q(7) on the atom + is defined as 


Q(r) = No(r) — N(7), (20) 


where No (7) is the number of electrons in the analogous orbitals of the neutral atom. 
The bond orders p,; have been calculated according to Coulson [11] and [13]: 


Drs = >, Ny Cr Cys (21) 
J 


The value of NV; and the summation over j are as in (19). When the summation over 
j is over orbitals occupied by z-electrons only, p,; is called mobile bond order. 


f) Energy values 


As mentioned above, the orbital energies obtained by solution of (5) and (8) are 
equal to the ionization energy for removal of an electron in the corresponding orbital. 

The total z-electron energy for the case when all occupied orbitals contain two 
electrons is 


EH = >, Hp"? + > a. (22) 


The summation should be carried out over all occupied space orbitals. 
For the special cases FAD- and DPNH-+*, where an electron is unpaired, the total 
a-electron energy is 


E,,= > HP + 1/2 Hse + ert Lees (23) 


The summation over j should be carried out over all doubly occupied space orbitals 
and e indicates the singly occupied orbital. 

We are interested in differences in total energy between DPNH + FAD (I), DPNH+ 
+FAD- (IT) and DPN+ + FADH, (III). The total energy, ZH, consists of the z-electron 
energy and the core energy, and in II, also the electrostatic energy, Hy, which is 
dependent on the relative orientation of the charged radicals. Assuming that the 
same interatomic distances occur in I and II, the core energy is approximately the 
same in I and II. The core energy of III is, however, different, partly due to 
changed interatomic distances and partly due to the formation of new o-bonds. 
It therefore seems possible to calculate only the energy difference between I and 
II with the present method. When calculating the energy difference 


H(1)~ £ (Il) = £,(1) — #,(1D) — By, (24) 


the electrostatic energy, H.,, is calculated on the assumption that the radicals are 
situated in two parallel planes, the O(9) of FAD- situated above N of the pyridine 


104 


ARKIV FOR FYSIK. Bd 17 nr 5 


ring of DPNH* as in Fig. 7, and that the charges are distributed according to Figs 2 
and 3. The distance between the planes is assumed to be 3.3 A. (Hausser and Murrell 
[14] assume that the corresponding distance between two similar free radicals is 
between 3 and 3.5 A.) 


Ill. Results 


As mentioned above, a calculation in the Hiickel approximation has first been 
carried out for each substance to get ¢ values for the first cycle of the SCF calcula- 
tion. The results of these Hiickel calculations are given together with the final 
results from the semi-empirical method according to Pariser and Parr [1] to get 
a comparison between the methods. This latter method is denoted the SCF approxi- 
mation. The Hiickel approximation for the free radicals DPNH* and FA D- is obtained 


Table 4. Orbital energies ¢ in a.u. obtained with the semiempirical method, denoted 
the SCF approximation, according to (5) and (8), and in the Hiickel approxima- 


tion. 
DPN+ DPNH FAD FADH, DPNH+ FAD- 
\CF Hiickel SCF Hiickel SCF Hiickel SCF Hiickel (free radical) (free radical) 
SCF SCF 
Occupied orbitals 
—0.661 —0.670 — 0.422 
—0.613 —0.695 =(.580 —0.626 = 0.3517 
—0.603 —0.635  —0.557 —0.628 (ivsyilis) tsa) — 0.846 — 0.322 
E93 —1.030 —0.515 -—0.602 -—0.541 -—0.591 —0.499 —0.556 — 0.745 — 0.276 
.700 —0.790 —0.487 -—0.540 -—0.483 —0.572 —0.429 —0.526 — 0.708 — 0.240 
.694 —0.780 —0.414 —0.489 —0.449 —0.520 —0.390 —0.461 =a Or —OokOl* 
Empty orbitals 
265 —0.655 —0.006 -—0.370 -—0.088 —0.422 —0.021 —0.374 
.246 —0.585 +0.021 -—0.308 —0.016 —0.349 +0.019 —0.337 
* Only an unpaired electron in this orbital. 
+ 
DPN DPNH DPNH 
i N. ae 
“e087? C008 +0.00C . C t0.00 +0.09C C+0.09 
+012C C+0.12 -0.04C C-0.04 +0.13C ue C+0.13 
Se Dae See 
+0.11 | | 
-0.21H, +0.05H,, 


Fig. 2. The number 


(free radical) 


s give the charges on the atoms according to (20), with the electrons considered 
explicitly. The values are obtained with the SCF method. 
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directly by removal of an electron from the highest occupied level of DPNH and by 
transferring an electron to the lowest empty level of FAD. 

The orbital energies ¢ are given in Table 4 for the levels occupied in the ground 
state and for the empty levels of lowest energy. — €; gives the ionization energy for 
removal of an electron occupying the 7:th orbital. 


+0.04 +0.19 
N 
a oie eae b — ae 
MS Rene ey ya N+0.17 
| +0.27 +0.21 
H 0.26 


Fig. 4. a) Differences in gross atomic population between DPNH and DPNH+, calculated with 
the SCF method. b) Differences in gross atomic population between FAD- and FAD, calculated 
with the SCF method. 
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Table 5. The gross atomic population N (r) according to (19) with the electrons 
considered explicitly. 


. DPNH+ 
4 : DPN+ DPNH. (free radical) 

N(r) on atom SCF Hiuickel SCF Hiiekel SCF Hiickel 
C(1) and C(5) 

respectively 1.08 1.14 1.00 0.90 0.91 0.89 
C(2) and C(4) 

respectively 0.88 0.7% 1.04 1.10 0.87 0.90 
C(3) 0.89 0.49 0.83 0.62 0.87 0.54 
N (6) 1.20 1.78 1.88 1.66 1.63 12 
ETS (a) = — 2 1.72 0.95 1.62 

FAD- 
FAD FADH, (free radical) 
SCF Hiickel SCF Hiickel SCF Hiickel 

N(1) BLE 1.46 1.23 1.51 120 1.46 
C(2) 0.67 0.50 0.82 0.68 0.69 0.64 
N (3) 1.11 1.41 1.20 1.45 1.28 1.47 
C(4) 0.86 0.50 0.84 0.69 1.08 0.82 
C (5) 0.91 0.68 1.10 1.06 0.94 0.81 
C(6) 0.88 0.66 0.94 0.93 0.93 0.74 
N (7) 1.92 1.74 1.94 1.87 1.95 1.79 
N(8) 1.09 1.29 1.95 1.90 1.36 1.47 
O(9) 1.39 1.76 1.96 1.92 1.58 1.80 


The distribution of the electrons considered explicitly, i.e. the gross atomic popula- 
tions of these electrons, is given in Table 5. The charges @ on the atoms are given in 
Figs 2 and 3. The differences in the atomic populations between DPNH and DPNH* 
and between FAD~ and FAD are given in Fig. 4. The mobile bond orders according 
to (21) are given in Figs 5 and 6. 

The difference in z-electron energy between DPNH+FAD (I) and DPNH*++ 
FAD- (II) is — 0.24 a.u. This value is obtained with the semi-empirical method accord- 


DPN DPNH DPNH- 
BS iS oo Nee oo Neg 
Cy j 
ee ae C C é 
0.63 [oss 0.35] joss a. [os 
C C C C c C 
Woe Or. <P 
[ase Joss 
H, H, 


(free radical) 


Fig. 5. The numbers along the bonds give the mobile bond orders according to (21), with the SCF 
method. 
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[ox [o2s 


C N 
ne Dew 
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| 


Nea go Nte op OH 
C C 


0.44 


0.82 


C N 
| 
H 
FAD. 
I 
N fee 8? N Xe ol 0 
Cc C 
ozs 0.42 
Cc N Fig. 6. The numbers along the bonds give the 
Nea 2 0< mobile bond orders according to (21), with the SCF 


method. 
(free radi cal) 


ing to Pariser and Parr [1], here denoted the SCF approximation. The electrostatic 
energy for DPNH++FAD- is —0.12 a.u., if the distance between the planes is 
3.3 A (for the configuration in space see Fig. 7). The total energy difference between 
(I) and (II) is therefore —0.24 a.u.=—6.5 eV for large distances between the 
charged radicals and —0.12 a.u.= —3.3 eV for the distance 3.3 A, on the assump- 
tion that the core energy is unchanged and that no charge transfer complex is 
formed. The formation of such a complex is discussed in IV. 


IV. Discussion 


From Table 4 it is seen that the ionization energies obtained with the Hiickel 
approximation are larger than the values obtained with the semi-empirical method 
according to Pariser and Parr, denoted the SCF approximation. The differences 
between the values are especially large for the empty orbitals. For the occupied 
orbitals the agreement between the methods is comparatively good, the differences 
varying between 0.03 and 0.12 a.u. In the following discussion we shall only consider 
the values from the method denoted the SCF approximation. 


108 


ARKIV FOR FYSIK. Bd 17 nr 5 


< C9 Ni, 


C 
H, 
Fig. 7. Assumed orientation in two parallel planes 


of the free radicals DPNH+ and FAD-, when the 
electrostatic energy is calculated. 


It is seen that the ionization energies of an electron from DPN+ and from FAD 
are very large (0.694 and 0.449 a.u. respectively), for DPNH and FADH, the corre- 
sponding values being lower, i.e. 0.414 and 0.390 a.u. respectively. The orbital 
energies of the lowest empty levels of DPN+ and FAD are —0.265 and —0.088 a.u. 
respectively, the corresponding values for DPNH and FADH, being only —0.006 
and —0Q.021 a.u. respectively. These results correspond well with the facts that 
DPNH and FADH, are easily oxidized, i.e. they are good electron donors, and that 
DPN* and FAD are easily reduced, i.e. they are good electron acceptors. 

Pullman and Pullman [4] have only given K values for the orbital energies, where 
é=a+Kf. According to a private communication they estimate —/ to be approxi- 
mately 20 Keal/mole and at most 40 Keal/mole. If a value within this interval is used 
for —f and our value —0.4242 a.u. for «, the e-values obtained with Pullman’s 
K-values for DPNH, FAD and FADH, agree well with our SCF values for the occu- 
pied orbitals of highest energy and with our Hiickel values for the empty orbitals 
of lowest energy.! This is reasonable, since Pullman and Pullman have used a Hiickel 
approximation. For DPN*+ however the ¢-values obtained from the K values of [4] 
differ considerably from our values, the value for the highest occupied level being 
—0.46 a.u. and for the lowest empty level —0.41 a.u. if — / = 20 Keal/mole is used. 
The poor agreement may be explained by the difficulties in choosing parameter 
values for N+ in DPN*. 

Those atomic charges which are of special interest are the charges on the N atom 
of the pyridine ring of the diphosphopyridine nucleotide, and of N(8) and O(9) of 
the flavine. (Index numbers in parenthesis refer to Fig. 1.) 


1 The discrepancies at most 0.03 a.u. 
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Considering the diphosphopyridine nucleotide, it is seen that the charge on N 
increases by 0.25 units when an electron is removed from DPNH, so that the charge 
on N in the free radical DPNH+ becomes +0.37. In DPN? it is further increased 
by 0.43 to a value of 0.80. Pullman and Pullman [4] have calculated the electron 
distribution in the free radical DPN, i.e. when an electron is transferred to DPN?*, 
and they have found the difference in charge on N between DPN* and DPN to be 
0.18 units. 

Considering the flavine, it is seen that the negative charge on N (8) in FAD increases 
by 0.27 units, when an electron is transferred to the lowest empty orbital, so that 
the charge on N(8) in FAD~ is —0.36. The negative charge on O(9) in FAD in- 
creases at the same time by 0.19, so that the charge on O(9) in FAD™~ is —0.58. 
Pullman and Pullman [4] have found the increase in negative charge on N(8) to be 
0.28 and on O(9) to be 0.02. 

In [3] an hypothesis is made for the manner in which a high-energy phosphate bond, 
due to the oxidation taking place in the DPN-FA D-part of the respiratory chain, can 
be formed. The hypothesis is based on the assumptions that a charge transfer complex 
can be formed between DPNH* and FAD-, and that phosphate can be bound to the 
O(9) of FAD-, if at the same time the electrons of O~ in the phosphate ion are dis- 
placed towards the V atom of DPNH* (XN is supposed to be positively charged) so 
that the o-bond between the P atom and the latter O atom is weakened. It is then 
a significant result of the calculation that the positive charge on V in DPNH and 
the negative charge on O(9) in FAD are increased by transfer of an electron from 
DPNH to FAD. However, it is seen from the results that the energy of DPNH*+ 
FAD- (II) is higher than that of DPNH+FAD (1). If the first step of oxidation of 
DPNH by FAD is an electron transfer in the supposed manner, the energy of (II) 
must be lowered e.g. by formation of a charge transfer complex. In [15] Mulliken has 
made an approximate calculation on a pair of molecules, viz. an electron donor, 
benzene (denoted B), and an electron acceptor, iodine (denoted A). In this case the 
energy difference between A~ + B+ and A+B is 8 eV, when the components are a 
large distance from each other, and between 4 and 5 eV when the distance between 
the components is 3.4 A. Mulliken maintains that his calculations give considerable 
support to the theory that the formation of a charge transfer complex causes the 
observed absorption at A 2900. Certainly this calculation refers to other molecules 
than the present ones, but since the energy differences are of the same order of 
magnitude as the energy differences found in the present calculations between I and 
II, it seems to be very probable that in the present case overlap between the states I 
and IT may cause a lowering of the energy by a complex formation. The present 
calculation thus gives considerable support to the hypothesis in [3]. 
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